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Abstract. We present new R filtered CCD observations 
of 5 faint and moderately faint asteroids carried out be- 
tween October, 1998 and January, 1999. The achieved ac- 
curacy is between 0.01-0.03 mag, depending mainly on the 
target brightness. The obtained sinodic periods and am- 
plitudes: 683 Lanzia - 4l'6±0'?2, 0.13 mag; 725 Amanda - 
> 3^*0, >0.40 mag; 852 Wladilena - 41^62 irt'Ol, 0.32 mag 
(December, 1998) and 0.27 mag (January, 1999); 1627 Ivar 

- 4^'80 ± Ol'Ol, 0.77 mag (December, 1998) and 0.92 mag 
(January, 1999). The Near Earth Object 1998 PG un- 
ambiguously showed doubly-periodic lightcurvc, suggest- 
ing the possibility of a relatively fast precession (Pi = ll'3, 
P2=5'?3). 

Collecting all data from the literature, we determined 
new models for 3 minor planets. The resulting spin vec- 
tors and triaxial ellipsoids have been calculated by an 
amplitude-method. Sidereal periods and senses of rota- 
tion were calculated for two asteroids (683 and 1627) by 
a modified epoch-method. The results are: 683 - Ap = 
15/195 ± 25°, /3p = 52 ± 15°, a/b = 1.15 ± 0.05, b/c = 
1.05 ±0.05, Psid=0^1964156±0^0000001, retrograde; 852 

- Ap = 30/210 ± 20°, /3p = 30 ± 10°, a/b = 2.3 ± 0.3, b/c = 
1.2 ± 0.2; 1627 - Ap = 145/325 ± 8°, /?p = 34 ± 6°, a/b = 
2.0 ± 0.1, b/c = 1.09 ± 0.05, P,.»d=0^1999154±0^0000003, 
retrograde. The obtained shape of 1627 is in good agree- 
ment with radar images by Ostro et al. (1990). 

Key words: solar system - minor planets 



1. Introduction 

Ground-based modelling of the shape and the rotation 
of the minor planets requires high precision and long- 
term photometric observations. With the advent of the 
CCD era it has become possible to study much fainter 
minor planets than previously. The photometric methods 
of modelling are based on multi-opposition lightcurves of 
full phase coverage obtained at very different longitudes 
(De Angelis 1993, Detal et al. 1994, Szabo et al. 1999). 
Another important aspect is to detect possible coUisional 
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effects, e.g. multiperiodic lightcurves due to binarity or 
precession, as they may yield insights into the recent solar 
system evolution. 

We started a long-term observational project ad- 
dressed to photometric monitoring of selected minor plan- 
ets. The observing programme consists of asteroids with 
available multi-opposition lightcurves enabling applica- 
tion of different photometric methods in order to model 
their shape and rotation. First results of this project have 
already been published in Sarneczky et al. (1999) and 
Szabo et al. (1999). The main aim of this paper is to 
present new CCD observations carried out between Oc- 
tober, 1998 and January, 1999 and models for 3 minor 
planets. Observations, their limitations and applied meth- 
ods are discussed in Sect. 2, while Sect. 3 deals with the 
detailed observational results. 



2. Observations and modelling methods 

We carried out Rc filtered CCD observations at 
Piszkesteto Station of Konkoly Observatory on ten 
nights from October, 1998 to January, 1999. The data 
were obtained using the 60/90/180 cm Schmidt-telescope 
equipped with a Photometries AT200 CCD camera 
(1536x1024 KAF 1600 MCII coated CCD chip). The pro- 
jected sky area is 29'xl8' which corresponds to an angular 
resolution of I'.'l /pixel. 

The exposure times were limited by two factors: firstly, 
the asteroids were not allowed to move more than the 
FWHM of the stellar profiles (varying from night to night) 
and secondly, the signal-to-noise (SN) ratio had to be at 
least 10. This latter parameter was estimated by compar- 
ing the peak pixel values with the sky background during 
the observations. The journal of observations is summa- 
rized in Table 1. 

The image reduction was done with standard IRAF 
routines. The relatively high electronic noises and low an- 
gular resolution did not permit the use of psf-photometry 
and that is why a simple aperture photometry was per- 
formed with the IRAF task noao.digiphot.apphot.qphot. 
Unfortunately other filters were not available during the 
observing run and consequently we could obtain only in- 
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Table 1. The journal of observations, (r geocentric distance; A - heliocentric distance; A - ecliptic longitude; (3 - 
ecliptic latitude; a - solar phase angle; aspect data are referred to 2000.0) 



Date 


RA 


Decl. 


r(AU) 


A(AU) 


A 


13 


a 


683 Lanzia 


















1998 12 14/15 


00 12.78 


+19 


58.9 


3.25 


2.82 


20 


18 


17 


1998 12 16/17 


00 13.78 


+19 


49.5 


3.25 


2.84 


20 


18 


17 


725 Amanda 


















1999 01 26/27 


06 31.61 


+27 


13.5 


2.34 


1.44 


110 


24 


12 


852 Wladilena 


















1998 12 12/13 


11 40.37 


+27 


28.2 


2.98 


2.71 


163 


23 


19 


1998 12 14/15 


11 41.67 


+27 


31.9 


2.98 


2.68 


163 


23 


19 


1998 12 16/17 


11 42.89 


+27 


36.2 


2.98 


2.65 


163 


23 


19 


1999 01 24/25 


11 47.52 


+30 


52.9 


2.95 


2.18 


170 


19 


14 


1627 Ivar 


















1998 12 14/15 


05 03.41 


+10 


30.4 


2.22 


1.26 


76 


-12 


6 


1998 12 15/16 


05 02.00 


+10 


32.6 


2.23 


1.26 


76 


-12 


6 


1998 12 16/17 


05 00.62 


+10 


34.9 


2.23 


1.27 


76 


-12 


7 


1999 01 22/23 


04 30.95 


+13 


09.2 


2.35 


1.65 


85 


-13 


20 


1998 PG 


















1998 10 23/24 


23 47.69 


+09 


15.0 


1.23 


0.26 


2 


9 


25 


1998 10 26/27 


23 55.11 


+08 


26.9 


1.23 


0.27 


2 


7 


25 


1998 10 27/28 


23 57.63 


+08 


11.6 


1.23 


0.27 


2 


7 


26 



strumental differential R magnitudes in respect to closely 
separated comparison stars. The precision was estimated 
with the rms scatter of the comp.— check magnitudes (tipi- 
cally 0.01-0.03 mag). 

We have also investigated the possible colour effects 
in neglecting standard photometric transformations. We 
made an R filtered 60-seconds CCD image of open cluster 
M67 on December 14, 1998. This cluster contains a widely 
used sequence of photometric standard stars (Schild 1983). 
We determined the instrumental magnitude differences 
in respect to star No. 81 in Schild (1983), which vari- 
ous colour indices arc close to zero {{B ~ V) = —0.098, 
(V — R)c = —0.047 mag). The studied standards were 
stars No. 106, 108, 117, 124, 127, 128, 129, 130, 134 and 
135, following Schild's notation. We plotted the resulting 
differences (Ai?j„s - ^Rstd) vs. {B - V) and {V - R)c 
in Fig. 1. For a wide colour range they do not differ more 
than 0.1 mag, while the colour dependence is quite weak. 
Therefore, the obtained instrumental i?-amplitudes of mi- 
nor planet lightcurves are very close to the standard ones, 
allowing reliable comparison with other measurements. 

The presented magnitudes throughout the paper are 
based on magnitudes of the comparison stars taken from 
the Guide Star Catalogue (CSC) (Table 2). Therefore, 
their absolute values are fairly uncertain (at level of ±0.2- Fig. 1. The colour dependence of instrumental minus 
0.3 mag). Fortunately it does not affect the other pho- standard magnitude differences for selected photometric 
tometric parameters needed in the minor planet studies, standard stars in M67. 
such as the amplitude, time of extrema, or photometric 
period. The final step in the data reduction was the cor- 
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Table 2. The comparison stars. The typical uncertainty 
in the magnitude values is as large as ±0.2-0.3 mag. 



Date 


Comp. 


m(GSC) 


683 Lanzia 










1998 12 14 


CSC 


1182 337 


15. 


.3 


1998 12 16 


CSC 


1182 85 


14. 


.4 


725 Amanda 










1999 01 26 


CSC 


1887 1325 


12. 


.3 


852 Wladilena 










1998 12 12 


CSC 


1984 2286 


12. 


.7 


1998 12 14 


CSC 


1984 2516 


12. 


.0 


1998 12 16 


CSC 


1984 2496 


13. 


.8 


1999 01 24 


CSC 


2524 1778 


12. 


.6 


1627 Ivar 










1998 12 12 


CSC 


702 759 


12. 


.6 


1998 12 14 


CSC 


689 1331 


12. 


.8 


1998 12 16 


CSC 


689 2101 


12. 


.6 


1999 01 22 


CSC 


681 519 


13. 


.7 


1998 PG 










1998 10 23 


CSC 


1170 1119 


14. 


.2 


1998 10 24 


CSC 


1171 632 


14. 


.3 


1998 10 26 


CSC 


1171 1424 


14. 


.5 



rection for the light timeQ Composite diagrams were cal- 
culated using APCll by Jokiel (1990) and are also light 
time corrected. Times of zero phase are included in the 
individual remarks. 

Two methods were applied for modelling. The first 
is the well-known amplitude- method described, e.g., by 
Magnusson (1989) and Michalowski (1993). For this the 
amplitude information is used to determine the spin vector 
and the shape. An important point is that the observed 
A(a) amplitudes at solar phase a should be reduced to 
zero phase {A{0°)), if possible, by a simple linear trans- 
formation in form of A{a) = A{0°){1 + ma), m is a pa- 
rameter, which has to be determined individually and that 
can be difficult, or even impossible if there are insufhcient 
observations (Zappala et al. 1990). 

The other possibility is to examine the times of light 
extrema ("epoch-methods", "E-methods"). In this paper 
a modified version was used, which gives the sense of the 
rotation unambiguously. The pole coordinates can be also 
estimated independently. Further details can be found in 
Szabo et al. (1999) and Szabo et al. (in prep.), here we 
give only a brief description. 

The initial idea is that the prograde and retrograde ro- 
tation can be distinguished by following the virtual shifts 
of moments of light extrema (e.g. times of minima). From 
a geocentric point of view, a full revolution around the 
Earth causes one extra rotational cycle to be added (ret- 
rograde rotation) or subtracted (prograde rotation) to the 
observed number of rotational cycles during that period. 
The virtual shifts increase or decrease monotonically and 



their cumulative change is exactly one period over one rev- 
olution. Therefore, plotting the observed minus calculated 
(O— C) times of minima versus the geocentric longitude, 
we get a monotone function ascending or descending by 
the value of the period. The definition of the observed 
O— C is as follows: 



0-C = r„,i„ -{Eo + N-P, 



sid . 



, AT, 



(1) 



where Tmin means the observed time of minimum, Eq is 
the epoch, while the N integer number denotes the cycles 
(e.g. the number of rotation) of Psid period between the 
observed extremum and the epoch. AT means the time 
interval between Tmin and Eq, and () denotes fractional 
part. The theoretical O— C curve depends on the pole co- 
ordinates: 



O-C 



0-C 



id 



1 tan(A - 

— arctan — - — 

TT sm(pp 



Ap) 



B) 



(2) 



^ Individual data are available upon request from the second 
author (szgyOneptun . physx . u-szeged . hu) 



where A and B denote geocentric longitude and latitude; 
Ap and /3p are the pole coordinates. 

The main difference between the classical E-methods 
and this O— C method is that time dependence is trans- 
formed into the geocentric longitude domain. Because of 
the system's basic symmetries, the O— C diagrams are 
calculated for a half revolution and with the half side- 
real period. The fitting procedure consists of altering P^id 
until the observed times of minima do not give a mono- 
tone O— C diagram showing an increase or decrease of 
exactly 1. Fitting a theoretical curve (Eq. 2) to the ob- 
served points, the pole coordinates can be also estimated. 

3. Discussion 

683 Lanzia 

This minor planet was discovered by M. Wolf in Heidel- 
berg, on July 23, 1909. It was observed in the 1979, 1982, 
1983-1984, 1987 oppositions (Carlsson & Lagerkvist 1981, 
WeidenschiUing et al. 1990). Carlsson & Lagerkvist (1981) 
determined a rotation period of 41*322 and an amplitude of 
0.14 mag. On the other hand, WeidenschiUing et al. (1990) 
measured a period of 4^37 with an amplitude of 0.12 mag. 

Our observations in 1998 suggest a period of 4l'6±0l'2 
with an amplitude of 0.13±0.01. Composite diagrams cal- 
culated with previously published periods between 4^" 3-4^^4 
have much larger scatter. The light-time corrected com- 
posite diagram is presented in Fig. 2. The zero phase is 
JD 2451162.3169. 

Based on earlier data (see Table 3), a new model has 
been determined with amplitude method. The observed 
amplitudes vs. ecliptic longitudes with the fit are plot- 
ted in Fig. 3. The resulting triaxial ellipsoid has the fol- 
lowing parameters: a/b=1.15±0.07, b/c=1.05±0.05, while 
the spin vector's coordinates are Ap = 15/195 ± 25°, 
/3p = 52 ± 15°, respectively. We could not reduce the ob- 
served amplitudes to zero solar phase, since the actual 



4 



Kiss et al.: CCD photometry and new models of 5 minor planets 



683 Lanzia 
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Fig. 2. The composite R lightcurve of 683 (symbols: solid 
squares - December 14; open circles - December 16). 
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Fig. 3. The observed amplitudes vs. longitudes with the 
determined fit for 683 Lanzia 



value of m parameter (e.g. Zappala et al. 1990) could not 

be estimated by the data sequence or asteroid classifica- 
tion. Also we have to note that a mixture of V and R am- 
plitudes was used, thus the model should be considered 
as an approximate one. The O— C model has also been 
determined (Fig. 4). For reducing the errors, lightcurves 
obtained between October, 1983 and February, 1984, were 
composed and one time of minimum was determined from 
this composite lightcurve. The resulting sidereal period is 
Pgirf=0di964156±0^0000001 with retrograde rotation. 



725 Amanda 

It was discovered by J. Palisa in Vienna, on October 21, 
1911. To our knowledge, the only one photometry of 725 
in the literature is that of Di Martino et al. (1994) carried 
out in 1985. They determined a sinodic period of 3^749 

associated with a full variation of 0.3 mag. Our observa- 
tions do not exclude that period, as they suggest a possible 
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Fig. 4. The observed O— C values fitted with the model 
for 683 Lanzia 



Table 3. Published photometries of 683 Lanzia 



Date 


A 


P 


Q 


A 




rcf. 


1979 03 19,20 


182 


-27 


9 


0^12 


43963.452 


(1) 


1982 12 16 


49 


9 


11 


0.14 


45319.591 


(2) 


1983 10 12,13 


130 


-9 


18 


0.15 


45650.871 


(2) 


1983 11 15 


137 


-13 


19 


0.15 


45650.871 


(2) 


1984 02 21 


129 


-23 


10 


0.16 


45650.871 


(2) 


1987 10 19 


16 


23 


7 


0.12 


47118.538 


(2) 


1998 12 14,16 


20 


18 


17 


0.13 


51162.275 


p.p. 



References: (1) - Carlsson & Lagerkvist 1981; (2) - Weiden- 
schilling et al. 1990 



value around 4 hours. Unfortunately the data cover only 
3 hours, thus we could not draw a firm conclusion. The 
observations were made under fairly unfavourable condi- 
tions, which is illustrated with the comp-check curve bear- 
ing a relatively high scatter (about ±0.03 mag). It is pre- 
sented together with the observed lightcurve in Fig. 5. 

852 Wladilena 

This asteroid was discovered by S. Belyavskij in Simeis, 
on April 2. 1916. Its earlier photometric observations were 
carried out in 1977, 1982 and 1993 (Tedesco 1979, Di Mar- 
tino & Cacciatori 1984, De Angelis & Mottola 1995). The 
observed light variation in 1998 had an amplitude of 0.32 
mag, while the period was 4l'62±0l'01. This is in very good 
agreement with results by De Angelis & Mottola (1995), 
who found a period value of 4l'613. The light time cor- 
rected composite diagram is presented in Fig. 6. The zero 
phase is at 2451160.5904. The lightcurve has remarkable 
asymmetries - the brighter maximum is rather sharp, its 
hump is exactly two times shorter than the other one. 
There are also small amplitude, short-period humps on 
the longer descending branch. These phenomena can be 
more or less identified in the previous measurements too. 
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Fig. 5. The R lightcurve of 725 on January 26, 1999. 



852 Wladilena 
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Fig. 7. The R hghtcurve of 852 on January 24, 1999 



852 Wladilena 



Table 4. Published photometries of 852 Wladilena 




0.2 0.4 0.6 0.8 

rotational phase 

Fig. 6. The composite R lightcurve of 852 (symbols: solid 
circles - December 12; dotted circles - December 14; solid 
squares - December 16). 



That is why we carried out a second observing run on 
January 24, 1999. We wanted to check the reality of these 
irregularities. The lightcurve revealed the same asymme- 
tries as those of observed one month earlier (Fig. 7). This 
may refer to a shape with sharp asymmetries, e.g. some- 
thing similar to a jagged tenpin. 

We have tried to determine a new model using the 
earlier data summarized in Table 4. Unfortunately, the 
measurements have such a distribution along the longi- 
tude that reliable modelling is difficult. This is shown in 
Fig. 8, where the observed amplitudes vs. ecliptic lon- 
gitudes are plotted with an approximate fit. The result- 
ing parameters are as follows: a/b=2.3±0.3, b/c=1.2±0.2, 
Ap = 30/210 ± 20°, /?p = 30 ± 10°. The pole coordinates 
are in considerable agreement with those of by De Angelis 
& Mottola (1995), who determined two possible solutions: 
(1) Ap = 53 ± 6°, /3p = 24 ± 20° and (2) Ap = 235 ± 6°, 



Date 


A 


P 


a 


A 


ref. 


1977 02 14 


139 


31 


10 


1^12 


(1) 


1982 10 18 


6 


-10 


10 


0.37 


(2) 


1993 11 8,10 


33 


-8 


3 


0.23 


(3) 


1998 12 12-16 


163 


23 


19 


0.32 


p.p. 


1999 01 24 


170 


19 


14 


0.27 


p.p. 



References: (1) - Tedesco 1979; (2) - Di Martino & Cacciatori 
1984; (3) De Angelis & Mottola 1995 



852 Wladilena 




21 ±20°. 



-50 50 100 150 200 250 300 350 400 

longitude 

Fig. 8. The observed amplitudes vs. longitudes with the 
determined fit for 852 Wladilena 



1627 Ivar 

This Earth-approaching asteroid was discovered by E. 
Hertzsprung in Johannesburg, on September 25, 1929. 
There are four photometric observations in the literature 
(Hahn et al. 1989, Velichko ct al. 1990, Hoffmann & Gcycr 
1990, Chernova et al. 1995) and one radar measurement 
by Ostro et al. (1990). The previously determined peri- 
ods scatter around 4'?8, thus our resulting 4^'80±0.01 is 
in perfect agreement with earlier results. The amplitude 
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Table 5. Published photometries of 1627 Ivar 
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0.2 0.4 0.6 0.8 

rotational phase 



1.0 



Fig. 9. The composite R lightcmve of 1627 (symbols: solid 
squares - December 14; dotted circles - December 15; 
crosses - December 16). 
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Fig. 10. The R lightcurve of 1627 on January 22, 1999 



changed significantly over a period of one month, as it was 
0.77 mag and 0.92 mag in December, 1998 and January, 
1999, respectively. The composite lightcurve is presented 
in Fig. 9, while the single lightcurve obtained in January 
is plotted in Fig. 10. 

A new amplitude model has been determined after col- 
lecting all available data (Table 5). The observed ampli- 
tudes were reduced to zero solar phase. First of all, the m 
parameter was derived from our measurements. The ob- 
served amplitudes in December, 1998 and in January, 1998 
were compared. As the longitudes differ by only 10°, and 
the difference between the corresponding phases is quite 
high (13°), the amplitude change can be mostly associated 
with the phase change. The result is m = 0.018. We have 
also corrected other amplitudes to zero solar phase and fit- 
ted the amplitude variations along the longitude. The cor- 
responding parameters are: a/b=2.0±0.1, b/c=1.09±0.05, 
Ap = 145/325 ± 8°, /3p = 34 ± 6°. The reduced amplitudes 
with the determined fit is presented in Fig. 11. The relia- 



Date 



A 



ref. 





1985 06 13 


317 


29 


48 




46226.750 


(1) 




1985 08 31 


15 


-21 


32 


0.55 


46258.703 


(1) 




1985 10 16 


4 


-23 


20 


0.63 


46287.184 


(1) 




1989 05 01-23 


203 


25 


20 


1.0 


47647.402 


(2) 




1989 06 15-23 


201 


21 


51 


1.12 


47647.402 


(2) 




1989 07 14-19 


213 


14 


60 


1.45 


47721.565 


(2) 




1990 05 11-14 


204 


25 


24 


1.08 


48029.439 


(3,4) 




1998 12 14,16 


76 


-12 


79 


0.77 


51162.295 


p.p. 


1999 01 26 


87 


-13 


18 


0.92 


51201.171 


p.p. 



References: (1) - Hahn et al. 1989; (2) - Chernova et al. 
(3) - Velichko et al. 1990; (4) Hoffmann & Geyer 1990 

1627 Ivar 



1995; 
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longitude 

Fig. 11. The reduced amplitudes vs. longitudes with the 
determined fit for 1627 Ivar 



bility of this model was tested by a direct comparison with 
radar images of Ostro et al. (1990). This is shown in Fig. 
12, where we used Fig. 5 taken from Ostro et al. (1990) 
with kind permission of the first author. The similarity is 
evident. 

The O— C method was used to determine the side- 
real period and the sense of the rotation. The results 
arc P^,d=0'?1999154±0^0000003, retrograde rotation with 
Ap = 143 ±8°, = —37 ±6° pole coordinates. The agree- 
ment between the poles obtained by different methods is 
very good. The sidereal period agrees well with results 
of Lupishko et al. (1986) - 0^19991, prograde -, but the 
senses are in contradiction. The fitted O— C diagram is 
presented in Fig. 13. 



1998 PG 

The Near Earth Object (NEO) 1998 PG was discovered by 
the LONEOS project in Flagstaff, on August 3, 1998. We 
observed about 80 days after the discovery, in October, 
1998. We found complex, strongly scattering lightcurves 
(two of them are shown in Figs. 14-15), which did not 
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Fig. 12. A pole-on view of the photometric model (left) 
and radar profile (right) of 1627 Ivar. The small ticks cor- 
respond to the uncertainties of the fit. 



1627 Ivar 




50 100 
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200 



Fig. 13. The observed O— C values fitted with the model 
for 1627 Ivar 



show any usual regularity. Therefore, we performed a 
conventional frequency analysis by calculating Discrete 
Fourier Transform (DFT) of the whole datasct (Fig. 16). 
Data obtained on October 27 are too noisy, thus we ex- 
cluded them from the period determination. 

The determined periods are 1^3 and 5^3, although 
these values have large uncertainties (about 10-15%). As- 
suming that the shorter period is due to rotation, we get 
a rotational period of 2^6. We note that our period val- 
ues do not contradict those obtained by P. Pravec and 
his collaborators, who found Prot=2^'517 and P2 ~ 7^*0 
(Pravec 1998, personal communication). The reason for 
doubly periodic lighcurve can be precession and/or bina- 
rity. The observed rate of multiperiodic lightcurves among 
NEOs is quite high (see, e.g., Pravec 1999), but the under- 
lying physical processes can only be identified with more 
detailed observations than we have on 1998 PG. There- 
fore, we conclude that we may have found evidence for 
precession in 1998 PG, but other explanations cannot be 
excluded. 

We summarize the resulting sinodic periods, ampli- 
tudes and models in Table 6. 
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Fig. 14. The observed and fitted R hghtcurves of 1998 PG 
on October 23, 1998 
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Fig. 15. The same as in Fig. 14 on October 26, 1998 
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